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ABSTRACT 


This investigation was undertaken to determine the feasibility of 
computing the friction velocity from available micrometeorological data. 
This could be of great assistance in finding the wind velocity profile by 
Deacon's equation and the eddy shearing stress in the lower portion of the 
friction layer. 

This investigation was carried out at the U. S. Naval Postgraduate 
School, Monterey, California, during the period February 1959 - May 1959, 
in partial fulfillment of the requirements for the degree of Master of 
Science in Meteorology. 

The author wishes to acknowledge the assistance and guidance given 
him by Professor F. L. Martin. He also wishes to express his appreciation 
to his wife for her patience and understanding during the many nights 


spent by him with the electronic computor. 
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1. Introduction 

The friction velocity, a stability parameter, and a term related 
to the Richardson number weré@ computed using micrometeorological data. 

The equations used fer the computation of these three parameters are 
derived. These equations are unique in that they combine the different 
theories ef two well-known investigators of surface layer turbulence inte 
an equation for the friction velocity and an equation for the wind profile 
in the surface layer. 

The micrometeorological data used in this investigation is discussed, 
as well as the instruments used in obtaining the data and the observa- 
tional errors. A theory is presented on errors of aerated temperatures 
close to the ground. 

After modifying the derived equations into a form suitable for using 
with an electronic digital computor, the so-called 'constants'"’ in these 
equations are discussed. The advantages and disadvantages of using an 
electrenic computor for this type of an davesctsation are then summarized. 

Various figures and graphs are used to illustrate the diurnal varia- 
tion as well as the possible variation with height as determined from 
this investigation. A comparison of these results with those of previous 


investigators is made, indicating close agreement in most cases. 





2. Derivations 

Deacon [2], in the derivation of his wind profile equation, as well 
as most other well-known investigators in the field of surface layer 
turbulence, considered only two virtual turbulence parameters - namely, the 


mixing length £ and the friction velocity u,, where 


Here © is the eddy stress and pis the air density. Lettau {10] introduces 
* 

a third parameter, the vertical mixing velocity w, that is quite important 

if nonadiabatic rather than neutral conditions are investigated. He also 


+ 
has a term which he calls the friction velocity u, where 


However, when the friction velocity is mentioned in this paper it will refer 
te the first definition. 

The surface layer is defined as that layer, immediately above the 
earth's surface, in which the eddy stress is constant with elevation through- 
eut the layer. Eqs. (1-9) are descriptive of this issech The detection. is 
that employed by Lettau in [10]. 


In a neutral surface layer the equations below are applicable. 


* * 
u Zz ‘a 


a a (1) 


‘Thewugheut this paper u will be understood to mean the wind velocity 
at some level averaged over some increment of time - say six minutes. 





Hes = k(z + a) k + von Karman's constant (2) 
Zz, * roughness parameter 


7. z ae Quy (3) 


ee (4) 


“s = x 
seg! a = tar (5) 


xo 
It should be noted that in Eqs. (1-5), wie Wes and TU. are constant 
with elevation in the neutral surface layer. 


In a non-neutral surface layer the next four equations are used. 


* * 
Ww 


u é (6) 
aos 3 (7) 
y= fF (8) 
(ie © (9) 


In [10], Lettau derives a relationship between “ and a which is of 
fundamental importance in the development of the remainder of this paper. 
Two basic assumptions are necessary in the derivation of Lettau's relation- 
ship: 

1. The ratio w/2 does not depend upon variations of thermal stratifi- 
cation, so that 
w 


e s a (10) 
“k fa 





Eq. (10) indicates that the vertical mixing length is changed in exactly 
the same ratio as the vertical mixing velocity if the thermal stratifi- 
cation varies. 

2. The second assumption is that the buoyant acceleration influences 
only the vertical parameters of turbulence. This assumption is necessary 
in order to fix absolutely the variation of mixing velocity and mixing 
length with thermal stratification. In an adiabatic atmosphere the term 
wal, can be considered to represent the vertical component of a turbulent 


acceleration. Thus we have the vertical acceleration equation: 





Ww n/n -Y : (11) 
x Ay a (pe 


where g = acceleration due to gravity, 36/dz = change in potent idl tempera- 
ture with respect to height, and Th «= mean absolute temperature averaged 
with respect to height in the layer considered. Eq. (11) indicates that 

the mean vertical turbulent acceleration in a non-adiabatic surface layer 


equals the adiabatic acceleration plus the buoyant acceleration. 


Upon multiplying Eq. (11) by 5 LAW 





OO 72 
2 ¢ F5r 
We. ky Tr W2 


% 2/k2 2/5 
but from Eq. (10), W/w 3 = £1¢ , 80 upon substituting 


2 08 p2 
we . 2. Loree 
fz hy Tn We 


multiplying through by yf? gives 





L ¥*,2 
Tm 
C6 pe 
or, was -/<; T5740 
k Ten Wy 
2 
g 28 
Let x = Baza GO 
mae 
mw? 
Therefore, Xo 4 —/ , : (13) 


Thus we finish with Lettau's derivation by putting 4h, in terms of x, 


By fe (14) 


a (1+X) 


The parameter x is important in the computations to be made in this paper. 
The next derivation was originally developed by Martin [12] combining 
Deacon's wind profile with Lettau's vertical mixing velocity. 


Starting with Deacon's Equation 


au = Uy (4 2e 22)" (15) 


where fp» as defined by Deacon, is a decreasing function of the Richardson 


number. Combining Eq. (7) and (8), one obtains 
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* 
* 
where uz {£ QU 
dz 
a + 
and w= AW 
Ay 
So that u, = Rp au 
Ae oz. 
Upon substituting into Eq. (15) 
LINN, Qu ae 
ou = NA Eo oz 2tZo5 
OZ Ren Dix. 
But from Ey. (14) & = i , so that after multiplying inside 
dy Gia) 


the radical by Aa fbr, 


A AIO A IE ET RT RL ETE 


JES Wa Quy -/{? 
Ou VOlsx) (+x). ae ZtZo 
a2 RZ Zo 
Dividing through by OU 
eh. 





Squaring out: 





-273 
Pole cee 0) Zee 
eyA Real) ( Sp 
But from Eq. (2) Ee = R({z te) , so that 
gu Rive (z32.) tee Tae 
OZ Roam + x) Zo 
Or. 
du ee We (252) ee 
oz R2,C+x)? Ze ze, 


Further combining of terms gives the equation in its final form. Namely, 


Ce ge zezey (16) 
OZ  R2,(itx)? Zo 


Thus Martin has eliminated u, from Deacon's Equation - a parameter that 


is difficult to determine with accuracy. 


Rewriting Eq. (16) in the form 





Ou = We = ie wee 2 = ae (17) 
PY 4 (i+x)* Zo ez. ( z) 








and comparing it with Eq. (15), one can easily see that an equation for 


the friction velocity has evolved: 


x i- 
ug oe We 24%) a (18) 
(1+ x)? Ho 


Using this equation to find u, one needs only to know the adiabatic mixing 
velocity, roughness parameter, vertical temperature gradient, and the 
Stability parameter. 

This completes the necessary derivations. The remainder of the paper 
deals with the solution of Eq. ( 2) for x, Eq. (16) for f7 , and Eq. (18) 


for u,, using available micrometeorological data. 





3. Data 

Immediately upon deciding to conduct this investigation into micro- 
meteorology a search was initiated to find suitable data. It soon became 
apparent that there was a scarcity of micrometeorological data of the 
caliber necessary to conduct individual research into surface layer turbu- 
lence. Several sets of available published data had to be discarded 
because one or another necessary parameter had not been observed, the period 
of observation was too short, or the heights of the observations were not 
at sultable locations. 

The best data that could be found within a reasonable time was that by 
Gerhardt [5].in 1949 and by Gerhardt [6] in 1950. The 1949 data that was 
of interest in this investigation consisted of wind observations at 12', 
41', and 91'; aerated temperatures at 3', 6', 12', 20', 35', 55', 80', and 
110'; non-aerated temperatures at 1", 6", 18", 3', 5', 9', and 15'. The 
1950 data included in addition to the above, wind at 3' and non-aerated 
temperatures at 23', 33', 45', 60', 78', and 99'. 

The information was gathered during runs of from 24 to 44 hours for 
different seasons of the year, with an average run being about 30 hours. 
Observations were taken every half hour during a run, with the final value 
for that time and height being an average of five readings taken over a 
six-minute interval. 

Jehn [9] and Gerhardt [7] give a detailed explanation of the instruments 
used and their accuracy, the topography of the micrometeorological site 
area, and the vegetation at the site and surrounding fields. The wind 
measurements were taken with Bendix-Friez Aerovanes to the nearest mph. 

The errors involved were less than 1 mph for wind speeds in the range 


10-50 mph, increasing to nearly 2 mph after 13 months of operation. For 


wind speeds less than 10 mph the errors were slightly greater, and readings 


9 





below 3 mph were unsatisfactory. 

The wind measurements, because of their inaccuracy and the spacing 
of instruments, were a severe handicap in this investigation. However, 
it was hoped that by analyzing a large number of observations, the errors 
would be somewhat minimized. 

The temperatures were obtained from calibrated rod-shaped thermistors 
(ceramic resistors), about 0.75 inches long and 0.03 inches in diameter. 
All thermistors were fitted to a single calibration curve giving an 
accuracy over the range of o°c to 50°C of bs 0.1°¢. The aerated thermistors 
were mounted in anti-radiation housings, being ventilated by a conven- 
tional motor and exhaust fan arrangement at the top of the housing. The 
non-aerated thermistors were coated with an anti-radiation paint and, in 
later runs, were shielded from both solar and terrestrial radiation without 
reducing the natural ventilation available. 

Because of all the precautions taken in calibration and shielding, 
ene could expect that there would be little or no difference between the 
aerated and non-aerated temperatures. However, during this investigation 
it soon became evident that there were large differences (in some cases 
nearly 1°c) between the two elements, and in several instances one gave an 
increase in temperature with height while the other gave a decrease over 
the same layer. In order to determine which would be best suited for this 
investigation, a thorough analysis of both sets of temperatures was con- 
ducted for two runs. It was found that at elevations greater than five 
feet the non-aerated temperatures were higher than the aerated both day and 
night, with the differences usually .2-.6°C. for all pairs of thermistors. 
For elevations lower than five feet, however, it was discovered that the 


non-aerated thermistors were cooler during the early afternoon and warmer 


during the early morning, with the differences increasing with decreasing 


10 





elevation. This at first seemed puzzling, since it is contrary to what 
one would intuitively expect -- at least for early afternoon. However, 
the author believes that this can be reasonably explained. 

The magnitude of the vertical temperature gradient decreases rapidly 
with height near the ground, especially during early morning and afternoon. 
Now consider the ventilation system drawing in an equal volume of air from 
both sheet and below the level z. Because of the difference in tempera- 
ture gradient the absolute value of the temperature difference between 
z and z - 42 will be greater than that between z and z +az. This dif- 
ference will be greatest during early morning and early afternoon, and 
thus the aerated temperature will be warmer than the true temperature at 
z during the afternoon and cooler in the morning. 

Because of this effect, it was decided that the non-aerated thermistors 
would give a more representative vertical gradient. However, of the four 
runs used in this paper (Runs 8, 10, 16, and 17) only on Run 17 was it 
possible to use the non-aerated temperatures since on Runs 8 and 10 they 
were only observed to 15' and on Run 16 the 3' element was broken. Even 
so, it is felt that greater accuracy in the temperature measurements are 
not warranted because of the large errors involved with the wind 
velocities. 

Gerhardt, in his final evaluation of the Micrometeorological Research 
Project (from which the data used in this investigation was obtained), 
stated: 

The importance of an accurate wind gradient, particularly in the 

lowest layers, cannot be overestimated; 12 measurement levels below 

50 feet would not be excessive. Great care should also be taken to 

obtain instruments with as nearly identical response as possible. 

Although thermistors have certain undesirable characteristics, it is 


believed that the installation as set up was satisfactory, and that 
accurate and representative air temperatures were obtained. 
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4. Computations 

In this section, the working equations using finite differences in 
place of derivatives are illustrated. The evaluation of the supposedly 
"non-variable" parameters is discussed, and the advantages and dis- 
advantages in using an electronic digital computer for this type of com- 
putation is pointed out. 

The particular working equation used to compute x was obtained by sub- 


stituting Eq. (2) into Eq. (12) and substituting finite differences in 


place of the partial derivative 00/92 . This results in the working 
equation 
X = | sk2(z+2z0)* | 2-8, (19) 
WA (Z27 2,) Tm 


The terms inside the bracket are constant for the layer during most 
of the run, although a new value of “. was selected every 12 hours. It 
is easily seen from this equation that x is zero for neutral conditions, 
greater than zero for unstable conditions, and less than zero for stable 
stratification. Values of AES were tabulated. 

The next variable to find was [? . The partial derivative 0u/o2 was 
first calculated from the wind data using finite differences. Then taking 
logarithms of both sides of Eq. (16) and solving for 1G resulted in the 


following equation. 


a 
fe = é I+ 4rk2. + by 4ot) 94 (20) 
Lyte Vi ere 
2+Z, 3+z, 
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The only variable for a given level is the numerator of the third term 
inside the brackets. 

The final variable to determine for this investigation was the 
friction velocity. The working equation was developed by taking logarithms 
of both sides of Eq. (18) and then taking e to the power of the left side 


and the right side, giving 


. J (21) 
Uy, = ett: — (br, BEF. 


This equation looks more difficult to solve than the original equation, 
but is actually more adaptable to programming on our electronic digital 
computer. 

Prior to solving Eqs. (19), (20), and (21) the "constant" parameters 
(von Karman's constant k, adiabatic mixing velocity wo and roughness 
parameter z) had to be determined. 

According to [8], _ and ZO are easily obtained by plotting the wind 
speeds against ln z for adiabatic conditions, and drawing a best-fitting 
straight line. The roughness parameter being the intercept on the ln z 
axis and the mixing velocity being the slope divided inte yon Karman's 


constant, or 


22 

WA = h(U2- UW) LPNS, a 
tn 2» 
2) 


The first problem to be solved was that of determining the time, and 
therefore the wind-sounding, for neutral conditions. All of the layers 
examined up to 110 feet were of depth not exceeding 10 meters. In layers 


of such depth, bounded by thermistors, a dry adiabatic lapse rate 
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corresponds to a vertical temperature difference lot £0.1°C. In such 
shallow layers, a neutral lapse may be regarded as essentially iso- 
thermal. Special emphasis in this regard was given to the temperature dif- 
ference in the layer 0-15 feet above the ground. Runs 16 and 17 corre- 
spond to stabilities nearly neutral for the entire period. 

On clear days, a neutral surface layer is considered [8] to occur 
near sunrise or sunset. This was found to be the case, The actual time 
of isothermalcy usually occurred between half-hourly wind observations near 
sunrise. Hence determining appropriate wind values required the averaging 
of at least two consecutive sets of data centered about the time of iso- 
thermalcy. Actually it was decided to average wind speeds for an hour 
before and after the neutral case in order to reduce round off errors in 
the reported wind values, since these were given only to the nearest mph. 
These averaged wind speeds were then plotted against ln z and a straight 
line fitted. The straight line went through or nearly through all the 
plotted points. The values of z were consistent with those 12 heurs 
earlier or later. The adiabatic mixing velocity “ still varied from the 
morning to evening adiabatic case. This is, however, to be expected. The 
adiabatic mixing velocity thus obtained was used for six hours preceding 
and six hours after the time of neutral stratification. 

Suggestions have been made [16], [17] that k is a function of stability. 
However, most specialists in the field treat k as a constant independent of 
Stability. Values given for it in textbooks and in the literature range 
from 0.38 to 0.45, with most references giving k as 0.40. This is the 
value used in this thesis. The particular constant value used for k will 
not materially affect the results of correlating 7 with Cia 

It was decided to use the CRA 102-A general purpose automatic digital 


computer [19] rather than a slide rule or hand computer. This decision was 
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based on the large number of computations involved and the additional 
numbers of significant digits possible with the computer. 

There are two possible disadvantages in using such a computer, the 
first being that the data had to be converted from the decimal number 
System to the octal number system. If the data were whole numbers, such 
as the wind observations, this presented little difficulty because programs 
are available for the necessary conversion. However, when involved with 
mixed numbers, as in the temperature observations, it is necessary to con- 
vert the whole number part and the fractional part separately, and then 
add them together. It was possible, however, after a little experience to 
rapidly convert the whole numbers mentally, thus necessitating machine 
conversion of the fractions only. It was found that for a run of 24 hours 
with wind and temperature at four heights and the mean temperatures of the 
three layers (a total of 528 numbers) that the data could fhe converted and 
put on IBM cards in less than one hour. 

The second disadvantage was the writing of the programs for the com- 
puter. It is estimated that a total of 20 hours was spent writing the 
necessary programs (Ri, du, (le fp » u,, and a correlation coefficient 
program), and another to now spent testing these programs on the com- 
puter to make certain they were working properly. 

Once this preliminary work was accomplished the remaining work was 
relatively simple, since the programs were good for any run once a few con 
Stants were changed. An example of the time saved will be illustrated. 
Five friction velocities were computed by hand using logarithms and a slide 
rule, taking nearly one hour. Using the computer, the friction velocities 
for one 24-hour run and three layers (244 values) were computed, typed out 
in both octal and decimal, and IBN cards punched in 15 minutes: Assuming 


that the time for the hand calculations could be reduced by half, it would 


les 








still take over 24 hours or about 100 times as long as the machine. 

By having the computer punch IBM cards for all variables it was 
then an easy matter to determine correlation coefficients between the 
various parameters. It took the machine 36 seconds to compute the 
correlation coefficient between two parameters consisting of 48 observa- 
tions each. 

The computer can also plot the information taking a little over a 
minute to plot 48 values of two variables. This was especially useful 
because the relatively small range the various parameters had would make 
it very tedious to plot by hand. 

Once a program has been successfully tested, there is little chance 
of erroneous values. With hand calculations, arithmetical errors may 
occur at any time. 

One last advantage that seems worthy of mentioning is the fact that 
with all the information on cards, it would be easy for another person to 
use in an investigation along similar lines - such as in determining the 


heat flux or the coefficient of eddy viscosity. 
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5. Results 

A discussion of the results of this investigation and a comparison 
of these results with those of other investigators into surface layer 
turbulence will now be made. The numerical values of the parameters 
a A > and u, for each run may be found in tables in Appendices I, 
II, III, and IV. It will be seen in these tables that for Runs 16 and 17 
three layers were evaluated - 3 to 12 feet, 12 to 41 feet, and 41 to 
91 feet, Although Deacon's wind profile is usually considered applicable 
only to about 41 feet, a higher layer was analysed in order to determine 
the inconsistencies, if any, from those at lower elevations. Runs 8 and 
10 were used, even though it was only possible to get values for the 
middle layer (due to insufficient data), because of the much greater ranges 
of stability than for Runs 16 and 17. 

The diurnal variation of (l4x)", 2, and u,, as well as their 
correlation with one another is illustrated in Fig. 1, where the three 
parameters for Run 8 (12 to 41 feet) are plotted against time. Run 8 was 
taken on 27-28 September with clear skies and light to moderate winds 
prevailing for most of the run. [It is seen that ae generally follows 
the diurnal march of temperature, becoming greater than one after sunset 
and less than one after sunrise, whereas 7? has the opposite trend. For 
this particular run p was less than one for only a few hours just prior to 
midnight, indicating the layer was unstable for the majority of the run. 

The inverse relationship between @ and (as. can easily be recognized 
in Figs. 1, 2, and 3. The curves of Figs. 1 and 6 have been smoothed 


somewhat uSing a smoothing equation 


2 = 2ofE + -50 3 + 25/3 
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For a closer examination of the short period fluctuations of fp» see the 
Appendix. A negative correlation between 73 and Shey” is to be expected 
since x is directly related to the Richardson number 
2 
gu) 
x = 02 


aN 
OZ 


Ri 








The quantities f and tee in Figs. 2 and 3 are weighted means for the 
entire layer from 3 feet to 91 feet for each observation time of Run 17. 


The equation used to determine this mean being 
Vee | 29/3 oe 
ie a 88 ie Li a ey 


and a similar equation for aoe Fig. 3 is a scatter diagram of 2 and 
(l+x) , with the line fitted by eye. The scatter of the points is prob- 
ably mainly the result of the inaccuracy of the original wind measurements 
which were given only to the nearest mph. For Run 17 the linear 
correlation coefficient between 2 and Goer of -.54 was determined. This 
significant negative correlation is considered to indicate the soundness 

of the Lettau-Martin theory. 

In 1948, Sutton [18] deduced that the Richardson number must be 
independent of height. However, Deacon, in the same year, found from his 
observations and calculations that Ri varied almost linearly with height. 
Fig. 4 shows (14x)? for each layer plotted against time for Run 16. A 
decrease of Cee with height (and thus a decrease of Ri with height) is 
very apparent in this figure. The average value of CESS) was determined 


for each layer, using the usual equation to find a mean 
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The average values were then plotted against height, resulting in Fig. 5. 
The linear decrease of (140) 6 with elevation, for this run, substan- 
tiates Deacon's conclusions. 

In comparing the valies of 7 obtained in this investigation with those 
of previous investigators, they were found to fall within the same range. 
Calder [1] had values of 6 from .75 to 1.4, Deacon (1948) found that @ 
varied from .75 to 1.2, and Priestly and Sheppard [15] determined the 
range of (3 to be from .70 to 1.3. A minimum value of equal to .84 
occurred in this investigation during Run 10, with the maximum value of 


1.4 appearing during Run 8. 





ee 


Table of ranges of (@ as determined by the 
indicated investigators. 


Table 1 
re was found to vary considerably from one observation time to the 
next, with a maximum variation of around .2 not uncommon for Runs 8 and 
10. Longley [11] derived an equation. for 3 that made the usual assumption 
of zy, constant with changes in stability unnecessary. In this equation 
pP depends only on wind speeds at three levels if the heights are related 
by 2 = ole Jas Zar’, where r is any ratdo. His equation was 
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Using this equation, he computed values of @ from Deacon's data and found 
an average difference from Deacon's f of from 0.1 to 0.2 depending on the 
heights selected. He also determined (2 from other data for various times 
during both day and night. He found that for daylight hours the range of 
four values of taken during the same hour might lay between 0.1 and 0.2, 
and for inversion conditions the ranges may even be greater. His final 
conclusions were (1) a possible error of 0.1 or 0.2 does not permit dif- 
ferent types of stability to be clearly distinguished in many cases, and 
(2) it is impossible to calculate (3 with sufficient accuracy to warrant 
its use in any practical problems. Deacon, in answering Longley's article, 
agreed that (2 is difficult to measure accurately, and said, 

I have never been so optimistic as to suppose that in practical 

problems it would be profitable to attempt the routine measure- 

ments of 2 in the field. 

It should be remembered that Longley's method of measurement was 
based essentially on curve fitting and had no way of separating out the 
effects of instability. Longley's conclusions may no longer be tenable 
since in the present paper, use of (l+x) accounts for a considerable part 
of the variability of /3 and of u,. 

In perusing the available literature on surface layer turbulence, etc., 
no investigations into the possible variation of (3 with height were found. 
- is generally considered constant with height, at least in the lower 
40 feet. However, Deacon (1948) felt there was a possibility that f? 
might vary slightly with elevation, and that if it did, it should slightly 


increase. Fig. 6 indicates a possible variation of with height. 
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Figure 7. The variation of & with the logarithm of height for Run 16. 
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VE , from Run 16, was averaged for each layer using the same equation 


as that for (14K) a with the results as below. 


Veg , = 0.9702 

Z5 

avy 

Ee 1.0040 
y be) 


* 


A 1.0271 
a 


VS 


wn 


3) 


It was found that in plotting these values against ln z, as in Fig. 7, 
the points could be connected by a straight line. This suggests that 
increases with the logarithm of height in the friction layer. The average 
Gs for each layer was also determined for Run 17. Here it was found that 
the average /3 for the upper layer decreased slightly from that for the 
middle layer. However, a straight line connecting Weue for the middle 
and lower layer gave a similar slope as that found for Run 16. This could 
possibly illustrate the fact that Deacon's equation was not applicable to 
as great a height for Run 17 as for Run 16, due to a difference in thick- 
ness of the surface layer, or it could mean that the relationship of 
with height found from Run 16 was only accidental, 

Prior to this paper, there has been no equation for accurately determin- 
ing the friction velocity for non-neutral conditions without measurements of 
the surface drag or short period wind velocity fluctuations. Panofsky [14], 
by making several assumptions, derived an equation for u, which depends 
only on the wind speed at three elevations. However, upon comparing his 


results for u, with those obtained by measurements of the short period 


w 
velocity fluctuations, they were found to be less than satisfactory. 


Deacon [3], while making an investigation into the variation of the 
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eddy stress with height, found that values of u, at 75 feet were much 


Es 
more variable than those at five feet. He attributed this mainly to 
accelerations associated with large scale eddies. He also found that 

there seemed to be no systematic variation of u, vith height, although 
individual measurements showed large departures. 

Fig. 8 is the friction velocity from Run 16 for the three layers plotted 
against time, as determined by Eq. (18). Here it is seen that the fluctua- 
tions of u, for the highest layer are generally of greater amplitude than 
for the lower two layers. Although this run seems to indicate a possible 
increase of u, with height, Run 17 indicated a possible decrease with height. 


* 


However, there does not seem to be a systematic increase of u, between the 


layer 3-12 feet as compared ta 12-41 feet. Any possible conclusions that 


can be made here regarding the variation of u, in the layer 41-91 feet are 


* 
tentative, inasmuch as the Deacon profile may not be applicable to these 


elevations. However there is no clear-cut increase of u, between the layer 


12-41 feet and 41-91 feet. 
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6. Summary and Conclusions 
Martin's combination of Lettau's vertical nixing velocity with 


Deacon's wind profile equation gives values of /3 and u, which compare 


* 

favorably with those of other investigators. However, the equations con- 

taining the term ae will not apply for x < -l, which may have the 

significance of placing an upper limit on the depth of the surface layer. 

During this investigation x was never found to be less than minus one. 
The equatien for the friction velocity, Eq. (18), may become an 


accurate method for determining u, if it becomes possible to evaluate 


* 


with sufficient accuracy. Additional work is needed to compare u as 


*? 
determined from Eq. (18), with that obtained from surface drag measurements 
or short period wind velocity fluctuations, before an accurate evaluation 

of the equation can be made. 

By comparing a large number of cases it was found possible to use data 
that was otherwise undesirable because of the inaccuracies of the wind 
measurements. However, to determine the various parameters for any one 
observation, it is felt that wind measurements of much greater accuracy are 
needed. 

One of the most important results of this investigation is the determin- 
ing of a possible variation of 6 vith height and the correlation of @ with 


a). However, further work with a larger number of layers is needed to 


verify this result. 
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